The Ratio of W + N jets To Z°/Y + N jets Versus N As a Precision Test of the 

Standard Model 
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We suggest replacing measurements of the individual cross-sections for the production of W+ N 
jets and Z+ N jets in searches for new high-energy phenomena at hadron colliders by the precision 
measurement of the ratios (W+0 jet)/(Z+0 jet), (W+f jet)/(Z+l jet), (W+2 jets)/(Z+2 jets),... 
(W+N jets)/(Z+N jets), with N as large as 6 (the number of jets in ttH). These ratios can 
also be formed for the case where one or more of the jets is tagged as a b or c quark. Existing 
measurements of the individual cross sections for W — + eu+ N jets at the Tevatron have systematic 
uncertainties that grow rapidly with N, being dominated by uncertainties in the identification of jets 
and the jet energy scale. These systematics, and also those associated with the luminosity, parton 
distribution functions (PDF's), detector acceptance and efficiencies, and systematics of jet finding 
and b-tagging, are expected to substantially cancel in calculating the ratio of W to Z production in 
each N-jet channel, allowing a greater sensitivity to new contributions in these channels in Run II 
at the Tevatron and at the LHC. 
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INTRODUCTION 

The signatures of the leptonic decays of the heavy 
gauge bosons W or Z° accompanied by jets, W — > lv+ 
jets and Z° — ► £ + £~ + jets, are among the preeminent- 
eminent search channels in very high energy particle col- 
lisions for 'new physics', i.e. interactions or particles that 
are not part of the Standard Model 0,Qi0>Q|- Many ex- 
tensions of the SM predict new particles which have elec- 
troweak (EWK) couplings and decay into the SM gauge 
bosons W, Z° , and 7, accompanied by jets. For example, 
searches have been made in the W or Z° +jets chan- 
nels for supersymmetric particles |5l la, technicolored 
hadrons 0, heavy W and Z' bosons yaM [HI that might 
arise in extended gauge groups or from excitations in ex- 
tra spatial dimensions, c harg ed Higgs bosons |lllll2l |. and 
leptoquarks |I3J, [1J, [lj, [16| , among others. More gener- 
ally, any production of new heavy particles with quantum 
numbers conserved by the strong interaction and EWK 
couplings is likely to contribute to signatures with one 
or more EWK gauge bosons; additional jets will always 
be present at some level from initial-state radiation, and 
may also be created in cascade decays of new heavy par- 
ticles or from the decay of associated heavy particles. 

Within the SM, the top quark was discovered and its 
mass measured in the W+3/4 jets channel in which at 
least one jet was identified as a b-quark [TtI flit HH Eo f . 
The W+ 2 jets channel with b-quark identification has 
been used to search for the Higgs boson [2l| and for single 
top (tb) production in the W + bb signature [22^. Asso- 
ciated Higgs production via ttH is expected to produce 
W+6 jets, of which 4 are b-quarks; associated W and Z 
production via ttW or tiZ will also produce W+6 jets, 
of which two will be b-quarks. 

Precise measurements of the W + N jets _23J and 
Z° j<y* -|_ TV jets [iij channels, where N is the number of 
jets, for values of N between and at least 6, including 



the cases where pairs of the jets are either bb or cc, would 
thus provide a broad search in a number of possible sig- 
natures of physics beyond the SM. The importance of 
calculating the cross sections for these channels has l ong 
been recognized [S3, EH E3, H^, H3, H3, H3, US, EM 
I35I l3rl l37l 13^: the development of sophisticated Monte 
Carlo programs capable of handling more particles in the 
final state at leading order (LO), or in some cases, next- 
to-leading order (NLO), now enables us to contemplate 
much more precise tests in the upcoming Tevatron Run 
II and at the LHC. 

However, direct measurements of the production cross 
sections of W + N jets or Z° /j* + N jets signatures suf- 
fer from inherent theoretical and experimental uncertain- 
ties associated with the definition and measurement (and 
hence counting) of jets. Among the dominant experimen- 
tal uncertainties are the energy response of the detector 
to a jet ('energy scale'), additional energy contributions 
from the underlying event (that part of the pp collision 
not directly involved in the hard parton-parton collision 
that produces the W or Z), backgrounds from misiden- 
tified non-electroweak events, and jet acceptance. These 
effects and others can change the number of jets mea- 
sured in a given event. Uncertainties in the theoretical 
SM predictions are dominated by the choice of Q 2 scale, 
the parton distribution function (PDF), initial/final state 
radiation (ISR/FSR), and the non-perturbative evolu- 
tion of partons into on-shell particles that would then 
be detected. All of these effects combined mean that the 
measurement of a specific exclusive N-jet channel such as 
W+ 4 jets will be completely dominated by systematic 
uncertainties at the Tevatron in Run II and at the LHC. 

In this note we use the Monte Carlo programs Mad- 
Graph [H |H and MCFM [U |H [H to explore using 
the measured ratios of W + N jets to Z°/y* + N jets at 
each value of N to provide a much more precise test of 
the SM than can be made by measuring the cross sec- 



tions themselves |44j. The W bosons are assumed to 
be identified by the leptonic decay W + — > e + v, and 
the Z°/7* intermediate state by Z°/j* — > e + e~. In 
most of the above models of new physics the production 
of new particles decaying into W and Z° /j* + jet final 
states would change the ratio from its SM prediction. 
The uncertainties listed above, except the misidentifi- 
cation backgrounds, are expected to cancel to a large 
degree, and the backgrounds can themselves be made 
to partially cancel by deriving the iv and £ + £~ event 
samples from a common inclusive high-pr lepton sam- 
ple |45j. We use data from the CDF |46j collaboration 
from Run I at the Fermilab Tevatron to estimate the sen- 
sitivity to contributions from non-SM processes using the 
a(W + Njets)/a(Z°/~f* + iVjets) ratio method. 

EXPERIMENTAL UNCERTAINTIES IN 

a(W + JVjets) /cr(Z°/j* + TVjets) 

The CDF Collaboration has published comprehen- 
sive studies of inclusive [Jg] W — > ev + N jets and 
Z°/j* — > e + e~ + N jets production in pp collisions at 
^/s = 1.8 TeV |46J. The D0 collaboration has measured 
the ratio of cross sections (W+ 1 jet)/(W+0 jet) 0i as 
the D0 measurements are less extensive in the number of 
jets (N) and do not include measurements of Z° + jets, 
we focus here on the CDF measurements. 

The CDF W selection required an electron with E t > 
20 GeV and \r)\ < 1.0, and missing transverse energy |45j 
]firp > 25 GeV. The Z selection required one electron 
satisfying the same charged lepton requirements, and a 
second electron with E T > 20 GeV for \r]\ < 1.0, E T > 
15 GeV for 1.1 < |ry| < 2.4, and E T > 10 GeV for 
2.4 < < 3.6. Jet identification |49| was made with 
a cone size in i]-(f> space of AR — 0.4, a threshold of 
E T > 15 GeV, and an 77 range of |?7| < 2.5. Multiple 
jets were required to be separated from each other in 77- 
4> space by a distance AR > 0.52; the requirement that 
the electron be 'isolated' from other clusters of energy in 
the calorimeter also corresponds to requiring AR > 0.52 
between the electron and each jet 

The individual (exclusive) cross sections extracted 
from the inclusive cross sections measured by CDF for 
W + N jets and Z°/ 7 * +N jets versus the number of jets, 
N, are displayed in Tableland Figure^ after being mod- 
ified for comparison with MadGraph's VF + predictions by 
dividing the CDF cross sections for W + + W~ by two. 
The uncertainties have been calculated in two ways: as- 
suming no correlations (giving an upper bound for the 
uncertainty) and assuming complete correlation (giving 
a lower bound). The uncorrelated uncertainties at each 
value of N have been calculated by subtracting the un- 
certainties of higher values of N in quadrature, and are 
reported first in the table. This overestimates the un- 
certainties, but as the (N+l)th channel is typically only 



20% of the Nth channel the overestimate is not large. The 
correlated uncertainties at each value of N have been cal- 
culated by subtracting the uncertainties of higher values 
of N; these uncertainties are reported second in the table. 

The estimated CDF systematic uncertainties are bro- 
ken down according to the source of each uncertainty in 
Table |H] versus the inclusive number of jets. One can 
see that in general the quoted systematic uncertainties 
grow rapidly with N, as described in detail in Ref. 0] 
This is due to the difficulties of counting jets given the 
rapidly falling spectrum in Et and the uncertainties in 
measuring the energy of a jet, and, to a lesser extent, 
uncertainties in the position of the jet with respect to 
the limit in rj in the jet selection. In addition, energy de- 
posited in the calorimeter from the fragments of the colli- 
sion not directly produced by the 'hard' interaction that 
produced the boson, called the 'underlying event', con- 
tribute to the total energy measured in the jet cone, and 
can promote jets from below threshold to over threshold, 
changing the number of jets in the event. Similarly, mul- 
tiple interactions from separate pp collisions in the same 
bunch crossing |50j can contribute energy in the jet cone. 
There are smaller contributions from uncertainties in the 
acceptance for the leptons, for the 'obliteration' of a lep- 
ton by a jet (if a jet lands close to a lepton the lepton 
can fail the identification criteria), and uncertainties in 
the contribution from decays of the top quark. Lastly, 
the uncertainty due to backgrounds from processes other 
than vector boson production ('QCD background') grows 
with the number of jets. 

The largest uncertainty is from the jet energy scale. 
This uncertainty will cancel in the production of W+jets 
and Z° /7* + jets events to the extent that the spec- 
tra in Et, the distribution in 77, and the composition 
(e.g. quark versus gluon) of the jets in the two processes 
are the same [Hi]]. Figure [21 shows the spectra in 77 and 
Et generated with the MadGraph Monte Carlo program 
|40j at LO. Using the difference of the ratio of the fitted 
slopes of the Et distributions for W and Z production in 
Figure times a typical uncertainty in the Et scale of 
20% |46J at 20 GeV gives an estimate of the uncertainty 
in the ratio of 2%. The effect of the finite acceptance in 77 
for jets depends on the difference in the distributions in 77 
of jets in W or Z production; taking the difference shown 
in Figure times the estimated variation with rapidity 
in jet response gives an estimate of the uncertainty 
in the ratio of 1%. 

The second largest systematic uncertainty is from the 
effects of energy from the underlying event, which can 
'promote' a 3-jet event to being a 4-jet event, for exam- 
ple, by boosting a lower-energy jet above the jet-counting 
threshold in Et- We expect that the underlying events 
in W and in Z events should be very similar; studies of 
the underlying event in jet events |5J| predict that the 
contribution from the beam fragments, which could be 
different due to the different quark diagrams in W and in 
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TABLE I: The cross sections times branching ratios for W + + N jets and Z° fa* + N jets production (first two columns) 
extracted from the CDF measurements versus the number of jets, at y/s = 1.8 TeV. These are used to calculate the ratios of 
the W + + N jets to Z° fa* + N jets jet cross section times branching ratio (third column). Also shown are the (less robust) 
ratios of a(W + N]ets)/a(W + N + ljets) and a(Z°/fa + N jets) /a(Z° fa* + N + ljets) . The first uncertainty given is the 
uncorrelated uncertainty, while the second (in parentheses) is the correlated uncertainty. These uncertainties are derived as 
discussed in the text. 



N 


a W+ + Nj 


C Z+Nj 


0- w + +Nj /&Z+Nj 


O W+Nj 1 & W+N+lj 


Cz+Nj faz+N+lj 



1 

2 
3 
4 


1010 ± 54 (34) 
185 ± 25 (17) 
41.5 ± 8.7 (6.5) 
7.7 ± 2.3 (1.4) 
1.45± 1.15(1.00) 


185.8 ± 11.1 (6.7) 
35.5 ± 5.5 (3.9) 
7.8 ± 1.8 (1.34) 
1.6 ± 0.53(0.39) 
0.43 ± 0.17(0.17) 


5.43 ± 0.44(0.27) 
5.21 ± 1.06(0.75) 
5.32 ± 1.65(1.23) 
4.78 ± 2.14(1.46) 
3.37 ± 2.99(2.68) 


5.46 ± 0.78(0.53) 
4.46 ± 1.10(0.81) 
5.42 ± 1.98(1.30) 
5.28 ± 4.48(3.77) 


5.23 ± 0.87(0.60) 
4.55 ± 1.26(0.93) 
4.88 ± 1.97(1.46) 
3.72 ± 1.92(1.73) 



TABLE II: The systematic uncertainties in percent on the measured CDF inclusive W + N jet production cross sections for 
for N=l to N=4 (column 1) 0. The successive columns are the uncertainties in the cross sections due to uncertainties in: 
the calorimeter jet energy scale, the underlying event, QCD background to W identification, multiple pp interactions in a 
single event, the value of the maximum allowed \rj\ for jets to be counted, the W acceptance, 'obliteration' of an electron by 
the superposition of a jet, and contributions from the top quark. The larger error bar is quoted in the case of asymmetric 
uncertainties. 



N(Jets) 


Etj Scale 


Und Ev 


QCD Bkgd 


Mult Int 


VJ 


Acc 


Oblit 


Top 


> 1 


6.8% 


5.8% 


5.2% 


3.2% 


1.9% 


0.8% 


0.2% 


0.05% 


> 2 


11% 


9.8% 


5.4% 


7.2% 


3.7% 


1.0% 


0.3% 


0.3% 


> 3 


17% 


16% 


9.1% 


9.8% 


4.8% 


1.8% 


0.6% 


1.3% 


> 4 


23% 


21% 


15.8% 


14% 


5.5% 


3.5% 


1.3% 


0.5% 



Z production, are a small portion of the total. However, 
the energy per tower contributed by the underlying event, 
and hence the effect on 'promotion' of jets, can be directly 
measured in W + N jets and Z° fa* + N jets events. We 
consequently assume that this uncertainty will be negli- 
gible in the ratio. 

For higher (> 4) jet multiplicities QCD backgrounds 
become comparable to each of the above. The back- 
grounds in the Z° fa* channel are at the few percent 
level, and are measurable (and hence subtractable) by 
counting same-sign events. Previous studies of the back- 
grounds to inclusive W production by CDF [5_4J for selec- 
tion criteria similar to those used here have shown that 
the background is dominated by approximately equal 
contributions from leptons from heavy flavor production 
and misidentified hadrons. How well these can be mea- 
sured with the new Run II detectors is not yet known; 
the former can be measured with the silicon vertex de- 
tectors, and the latter can be measured by conventional 
background techniques. 

The next largest systematic uncertainty in the Run I 
CDF cross section, contributions from multiple pp inter- 
actions, should cancel identically in the ratio, as it is 
uncorrelated with the hard scattering. 

The remaining uncertainties due to acceptance, 'oblit- 
eration' of a lepton by a jet, and contributions from top 
decay, are at most at the few percent level 



MONTE CARLO PROGRAMS AND EVENT 
SELECTION CRITERIA 

We have explored the W to Z° ratios in pp collisions at 
the Tevatron energy of yfs = 1.96 TeV using the Monte 
Carlo programs MadGraph |4(J and MCFM >43:]. Sam- 
ples of W+(-> e+u) + N jets and Z°/fa (-> e+e") + 
N jets , for N up to 4, were produced at LO using Mad- 
Graph. MCFM was used to explore the ratios for up 
to 2 jets at NLO, and to understand the dependence of 
the ratios on the Q 2 scale and on the parton distribu- 
tion functions for up to 4 jets at LO. Jets are treated at 
the 'parton level' with kinematic selections applied to the 
4- vectors with no fragmentation or detector simulation. 

We consider only the production in first-order elec- 
troweak processes of the W+jets and Z+jets channels - 
i.e. production of boson + jets from the WW, WZ, and 
ZZ channels are excluded. We also exclude tt and tb pro- 
duction; the method proposed here should allow a more 
precise determination of the non-top W+jets production, 
the dominant background in the top channel, and hence 
should allow more precise measurements of the top quark 
mass and cross section. 

The selection criteria and strategy for W and 
Z° fa* events used in the Monte Carlo studies were de- 
veloped for the measurement of R, the ratio of inclusive 
cross sections R = a(W)/a(Z° fa* ) [Hif. To minimize 
systematic uncertainties in the ratio due to the trigger 
and lepton selection, both W and Z° fa* events are se- 
lected from a common sample of inclusive central high 
transverse momentum |45| (pt) leptons, with transverse 
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FIG. 1: a) The measured cross sections for the signatures 
W + (-^ e + v) + N jets and Z°/j* (— > e+e") + N jets versus 
the number of jets, N, in W + and Z° /j* production in pp 
collisions at */s = 1.8 TeV. The data are from the CDF |4g] 
collaboration and were originally reported as inclusive cross 
sections. In computing exclusive cross sections from these, the 
uncertainties have been calculated in two ways. The dotted 
error bars were calculated assuming no correlations (giving 
an upper bound for the uncertainty) and the solid error bars 
were calculated assuming complete correlation (giving a lower 
bound); b) the percent uncertainty in the W + + N jets and 
Z° /7* + N jets cross sections. The uncertainties shown are 
the lower bounds (corresponding to the solid error bars in plot 
a). The figure shows the rapid growth of the uncertainties 
with N, the number of jets. 



energy (Et) greater than 25 GeV and pseudo-rapidity 



(I77I) less than 1.0. The second lepton from the boson de- 
cay, either another charged lepton (from Z°/y* decay) or 
a neutrino (from W decay), is required to have Et > 25 
GeV; in the neutrino case this is implemented by requir- 
ing the missing transverse energy (E T ) to be greater than 
25 GeV. 

Jets are required to have Et > 15 GeV and to be 
within 1 7/| < 2.5. Our MC studies are at parton level, so 
that there are no considerations of cone size, energy scale, 
or acceptance corrections in the Monte Carlo numbers. 



THE PREDICTED RATIOS 

a(W + + TVjets) /a(Z°/y* + TVjets) 

The predicted LO ratios <j(W+ + TVjets) Aj(Z° /<y* + 
TVjets) are presented in Figure and in Table II I II To 
determine the final uncertainty on this ratio will take a 
full analysis of the Run II data set; in lieu of this we 
have made some simple assumptions to get an estimate 
of the sensitivity in cross-section for non-SM physics in 
each of the N-jet channels in Run II of the Tevatron. We 
assume that the jet energy response of the calorimeter 
will largely cancel for jets in Z°/y* events and W events 
as discussed below. We also assume the effects of the 
underlying event in Z° /<y* and W events will similarly 
cancel. These are the two largest contributors to the 
systematic uncertainties quoted in Ref. [lif . 



THE PREDICTED RATIOS 

a(W + Njets)/a(W + TV + ljets) AND 

o{Z°h* + 7Vjets)/a(Z°/ 7 * + N + ljets) 



While the ratios of cross sections a(W + Njets)/a(W + 
N + ljets) and cr(Z°/ 7 * + A^jcts) / 'a{Z° / 7 * + N + 
ljets) are much more difficult to measure precisely than 
the a(W+ + Njets)/a(Z°/y* + TVjets) ratios, we include 
the generator-level LO predictions for them here as they 
are often used in extrapolations in N to estimate back- 
grounds at large N, and also to measure the strong in- 
teraction coupling. These are reported in Table ITTT1 and 
shown in Figure 



THEORETICAL UNCERTAINTIES IN 

a(W+ + iVjets)/£r(Z /7* + Wjets) 

The two largest uncertainties in the predicted LO W + 
N jets and Z° /-f* + N jets cross sections are expected 
to be due to choice of Q 2 scale and parton distribution 
function (PDF). We investigate the dependence of the 
ratio on these two choices using MCFM. 
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FIG. 2: The plot on the left (a) shows the (normalized) jet r\ distributions for W + {— + e + u) + 1 jet (dashed) and 
Z°/7* (— + e + e~) + 1 jet (solid) events satisfying the selection criteria described in the text. The plot on the right (b) 
shows the corresponding jet Et distributions, log (dN/dEr) versus N. Both plots are predictions at LO using MadGraph 40]. 
Uncertainties in the ratios a(W + + Njets)/a(Z° /j* + iVjets) due to the uncertainty in the jet rapidity cut at r)=2.5 are 
estimated from the shapes in the left-hand plot, and those due to the uncertainty in the jet energy scale from the right-hand 
plot. 



TABLE III: MadGraph leading order predictions of cross sections times branching ratio for W + + N jets and Z° /j* + 
N jets production (first two columns) versus the number of jets, at y/s = 1.96 TeV, which are used to calculate the ratios of 
the W + + N jets to Z° /j* + N jets jet cross section times branching ratio (third column). Also shown are the (less robust) 
ratios of a(W + N 'jets) /a(W + N + ljets) and a(Z° /-y* + JVjets) /a(Z° j-f + N + ljets) (last two columns). 
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a W++Nj 


CZ+Nj 


a W+ + Ni 


/oz+Nj 


VW+Nj f&W+N+lj 


&Z+Nj f&Z+N+lj 





341.5 ± 0.5 


67.0 ± 0.2 


5.10 ± 


0.02 


8.11 ± 


0.06 


6.41 ± 


0.02 


1 


42.1 ± 0.3 


10.45 ± 0.01 


4.03 ± 


0.03 


5.08 ± 


0.05 


4.91 ± 


0.07 


2 


8.28 ± 0.05 


2.13 ± 0.03 


3.89 ± 


0.06 


4.93 ± 


0.07 


4.75 ± 


0.09 


3 


1.68 ± 0.02 


0.448 ± 0.006 


3.75 ± 


0.07 


4.71 ± 


0.09 


5.15 ± 


0.09 


4 


0.357 ± 0.005 


0.087 ± 0.001 


4.10 ± 


0.07 











Dependence on the Q 2 Scale 



Q 2 also are listed. 



= aW + + N jets evaluated at Q 2 , 
1 ' - aZ /-f* + N jets . The ra- 



The effect of the choice of Q 2 scale is expected to par- 
tially cancel in W + N jets and Z° + N jets produc- 
tion, as both proceed through a Drell-Yan-like process. 
We define W$(Q 
and, similarly, Zn(Q 2 ) 
tios of W and Z cross sections evaluated at = My 2 
and at Q 2 = M v 2 + P T , V 2 , W+{M 2 ) /W+(P 2 + M 2 ) 
and Zn(M 2 ) jZ^(P 2 + M 2 ) , are given in Table Hvl and 
shown in Figure El Changing the Q 2 scale affects the 
W cross sections by as much as 15% and affects the Z 
cross sections by as much as 12%. However, changing 
the Q 2 scale has much less effect on the predicted ratio 
a(W+ + Njets)/a(Z°/j* + Njets) ; which changes less 
than 2%, as shown in Figure El and in Table IIVI where 
the W/Z ratios evaluated at the two different values of 



Dependence on the Choice of Parton Distribution 
Function 



We have used the MCFM generator and a selection 
of parton distribution functions to investigate the de- 
pendence of the cross sections in the VF + +2jets and 
Z° +2jets channels. The cross sections calculated 
with the CTEQ3L, CTEQ4L 0, MRST98 [13, and 
MRSG95 |5i| distributions were compared to the results 
calculated with CTEQ5L, the default PDF. The results 
of the comparison are reported in Table IVl Figure0gives 
the ratio of cr(VF++2jets for PDF x) to a(W++2 jets for 
CTEQ5L), while Figure El shows the ratio of cr(W^++2 
jets) to cr(Z+2 jets) for a given PDF. For the four PDF's 
we chose, the changes in the W and Z cross sections them- 
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TABLE IV: MCFM predictions for the ratios W + + N jets and Z /"/* + N jets with different Q 2 scales (columns one and 
two), and the ratio R + with different Q 2 scales. R + is a(W + + N]ets)/a{Z° /'y* + iVjets) , and (^corresponds to Q 2 = My 2 , 
while Q 2 2 corresponds to Q 2 = M v 2 + Pt.v 2 - 



N jets 


aW + (Q^) / aW + (Q 2 ^) 


oZ" (Qi a ) / 


aZ" (Q 2 2 ) 


R + (Q^) / R + (Q 2 ^) 





0.999 ± .001 


1.000 ± 


.001 


0.999 ± .001 


1 


1.017 ± .003 


1.018 ± 


.002 


0.999 ± .002 
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1.075 ± .002 


1.066 ± 


.002 


1.009 ± .002 
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1.153 ± .004 


1.134 ± 


.002 


1.017 ± .004 
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FIG. 3: A comparison of CDF data (circles) with Mad- 
Graph (LO) predictions (triangles) for the ratio of produc- 
tion cross sections times leptonic branching ratios for the 
signature W + + N jets to the signature Z° /j* + N jets , 
a(W + + N jets) /a(Z°/j* + iVjets) , versus the number of 
jets, N, in W + and Z°/j* production at y/s = 1.8 TeV for 
the data and at yfs = 1.96 TeV for the predictions. The case 
where two of the jets are b-quark jets are also shown (inverted 
triangles) . The statistical uncertainties on the predictions are 
smaller than the symbols. 

selves range from +27% to -7% for the W's and +25% to 
-8% for the Z's, while the the range of the change in the 
ratio is from +1.5% to zero, a factor of ~ 20 smaller. 



exclusive W + N jets cross section (if instead the source 
feeds Z+ jets at the same crossection, the sensitivity will 
be larger by a factor of about 10 [59| . 

The estimates above of the systematic uncertainties 
on are on the order of several percent; an estimate 
based on the Run I CDF experience in measuring R is 
that 1% in that ratio may be achievable [Hij]. Statistical 
uncertainties would then be expected to dominate over 
systematics in Run II at the Tevatron for N greater than 
2. 

Making the assumptions that the new contributions 
are to the W cross section and not that of the Z, that 
the systematics on the ratio can be reduced with a much 
larger dataset from several percent to 1%, and that 
one uses only the electron modes of W and Z decays, 
we find the 1-sigma cross-section uncertainties on new 
physics shown in Table IVD The muon channel would be 
expected to double the statistics (and hence lower the 
uncertainties by \/2). 

Additional sensitivity can come from comparing ob- 
served with expected kinematic distributions or by look- 
ing for additional objects in the events. In particular, 
the production of a pair of b-quarks suppresses the cross 
section over that for light quark production by a large 
factor, in principle allowing a corresponding increase in 
sensitivity. Table IVIII shows the ratio of the QCD cross 
section for producing N jets including no b quarks to N 
jets including two b quarks, for W or Z production. How- 
ever standard model top production will provide a large 
background for non-standard model physics in these sig- 
natures. 



CONCLUSIONS 



SENSITIVITY TO NEW CONTRIBUTIONS TO 

THE a(W + N jets) /a (Z°/j* + iVjets) RATIO 

A non-Standard Model source of W+ jets or Z + jets 
would result in a measured deviation from the expected 
SM value of the R N = a(W + Njets)/a(Z° /j* + Njets) . 
Assuming that the contribution is to W+jets, we can 
(crudely) estimate the sensitivity to new physics in each 
of the W+ N-jct channels by multiplying the uncertainty 
on the ratio a(W + Njets)/a(Z° /j* + A^jets) by the 



The measurement of the production cross sections of 
the vector bosons and Z° in association with a num- 
ber (N) of jets is now a standard way of looking for the 
production of new particles or processes that are not de- 
scribed by the Standard Model. With the expected in- 
creased luminosities of Run II and the LHC, N can be 
quite large; processes such as associated production of a 
Higgs boson with a ti pair can produce W+ 6 jets (4 of 
which are b-quarks), for instance. Increasing the preci- 
sion of the comparison with Standard Model predictions 
is necessary, as there are truly difficult problems, both 
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FIG. 4: The ratio of cross sections times branching ratios, a(W + + iVjets) /a(W + + N + ljets) (left hand plot) and a(Z° /j* + 
iVjets) /a(Z° /j* + N + ljcts) (right hand plot) versus the number of jets, N, in W and Z production at \fs = 1.8 TeV for 
the data and at yfs = 1.96 TeV for the predictions. The data (circles) are from the CDF and D0 01 collaborations; the 
predictions are at leading order from MadGraph |40J. The MadGraph cross sections for when the jets are from gluons or light 
quarks are shown with triangles, while inverted triangles represent when two of the jets are from b-quarks. 



TABLE V: MCFM predictions for the ratios W++2jets and Z° +2jets with different PDF's. The PDF's that were compared 
are CTEQ3L, CTEQ4L, CTEQ5L, MRSG95, MRST98. Column one gives the ratio of cr(W + +2jet) at one of the PDF's to 
<r(W /+ +2jets) at CTEQ5L. Column two is the analogous Z° information. The third column is the ratio of -R + (2) at a specific 
PDF to R + (2) at CTEQ5L, where -R + (2)= a{W + + 2jets)/cr(Z°/7* + 2jets) . 



PDF X 


(aW£) / (aWt) 


(ffZx) / (<tZl) 


(R£) / ( r l) 


CTEQ5L 
CTEQ3L 
CTEQ4L 
MRSG95 
MRST98 


1.000 ± .000 
1.103 ± .002 
1.105 ± .002 
1.268 ± .002 
0.932 ± .001 


1.000 ± .000 
1.090 ± .002 
1.094 ± .002 
1.249 ± .002 
0.922 ± .001 


1.000 ± .000 
1.011 ± .003 
1.009 ± .003 
1.015 ± .003 
1.011 ± .002 



theoretical and experimental, in predicting the cross sec- 
tions for W + N jets and Z°/j* + N jets when N is 
large. 

Using the Monte Carlo generators MadGraph and 
MCFM at the parton level, and the published CDF data 
on W and Z + jets production, we have made initial 
estimates of the systematic limits on the precision that 
can be achieved in the measurement of the ratios of W 
to Z production, a(W + N]cts)/a(Z° /V + Njets) , a s 
a function of the number of observed jets, N. The re- 
sults indicate that the ratios are at least an order-of- 
magnitude less sensitive to experimental and statistical 
uncertainties than the individual cross sections. In par- 
ticular the ratios are more robust for large values of N, 
where the experimental uncertainties in the energy scale 
and contributions from the underlying event and multiple 
interactions lead to a rapid growth in the cross section 
uncertainty with N. 



With respect to the theoretical uncertainties, at N=2, 
for example, we find the uncertainty due to choice of 
Q 2 scale is a factor of ~ 8 smaller in the ratio a(W + 
iVjets)/fj(Z°/7* + iVjets) than in the individual W or 
Z cross sections. Similarly, the uncertainty due to the 
choice of PDF, largely driven by the u/d quark ratio, is 
smaller in the ratio by a factor of ~20. 

The experimental uncertainties in the cross sections, 
dominated by the uncertainty in the jet energy scale and 
contributions from the underlying event, are greatly di- 
minished by focusing on the ratio of W and Z° /-y* cross 
sections ratio than the cross sections themselves. In par- 
ticular the uncertainty due to uncertainties in the jet en- 
ergy scale, the contributions from the underlying event, 
multiple interactions in one event, etc. cancel to a high 
degree. We have here made estimates at the parton level; 
a full determination of these will require the new data and 
a full analysis; our initial estimates are that the ratios can 
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FIG. 5: The ratios a{W + +N jets at Q 2 = M v 2 ) to a(W+ + 
N jets at Q 2 = M v 2 + Pt,v 2 ) and a(Z° /j* + N jets at 
Q 2 = M v 2 ) to a{Z°h* +N jets at Q 2 = M v 2 + P T ,v 2 ). 
Changing the Q 2 scale significantly changes the cross sections, 
by up to approximately 15%. However the ratio of W to Z 
cross sections changes much less (see Table ITvl . 



FIG. 7: The ratios <r(W+2jets for PDF X) to <r(W+2jets 
for CTEQ5L) and <r(Z+2jets for PDF X) to a(Z+2jets for 
CTEQ5L). Changing the PDF affects the cross sections quite 
significantly, by up to approximately 25%. However the ratio 
of W to Z cross sections changes much less (see Table . 




TABLE VI: The cross section corresponding to a 1-sigma un- 
certainty in the W/Z ratio in 2 fb~ , and in 15 fb -1 . The 
bins up through N=4 use the cross sections of (i^l : the N=5 
and higher bins have been extrapolated using an exponential, 
with a factor of 4.8 for each successive jet. Note that the 
number of Z° —+ e + e~ events in each bin will be approxi- 
mately a factor of 10 smaller than the corresponding number 
of W events. Using the dimuon channel one can gain a factor 
of approximately v2 on these uncertainties. 



N jets 



Event and W Properties 


W/Z Ratio Method Reach 


N(Jets) 


aw 




Tnew 15 fb 1 





1896 pb 


20 pb (1.0%) 


20 pb (1.0%) 


1 


370 pb 


4.4 pb (1.2%) 


3.7 pb (1.0%) 


2 


83 pb 


1.5 pb (1.8%) 


0.9 pb (1.1%) 


3 


15 pb 


0.5 pb (3.5%) 


240 fb (1.6%) 


4 


3.1 pb 


230 fb (7.5%) 


95 fb (2.9%) 


5 


650 fb 


100 fb (16%) 


40 fb (6%) 


6 


140 fb 


50 fb (36%) 


18 fb (13%) 


7 


28 fb 


20 fb (78%) 


8 fb (29%) 


8 


6 fb 




4 fb (63%) 



FIG. 6: The ratios (R 

2 



M v 2 ) to {R+ at Q 2 = M v 2 + 



at Q 2 

Pt,v 2 ), where R+ = a(W+ + Njets) /a(Z° /-y* + JVjets) 

2 



Changing the Q z scale affects this ratio by ~ 2% while the 
individual cross sections change by more than 15%. 



be determined at the several percent level. This is a sig- 
nificant improvement over the present uncertainties on 
the cross sections themselves. 

It is a pleasure to acknowledge helpful conversations 



with Edward Boos, John Campbell, Jay Dittman, Lev 
Dudko, Keith Ellis, Michelangelo Mangano, Stephen 
Mrenna, Jon Rosner, and Tim Stelzer. Special thanks 
are due to Tim Stelzer and John Campbell for help 
with MadGraph and MCFM, respectively. This work 
was supported by the National Science Foundation, grant 
PHY02-01792. 
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(N jets) / (66+(N-2) jets) for W + or Z 1 ' / 7 * 


W + + bb + N jets /Z u /~{* +bb + N jets 


W++2i/W+bb+0y. 90.29 ± .96 
W + +3j/W+bb+ly. 54.72 ± .84 
W + +4j/W + &&+2j:37.58 ± 1.30 


Z+2j/Zfe6+0j 
Z+3j/Z6fe+lj 
Z+4j/Z6fe+2j 


58.84 ± .89 
33.94 ± .69 
22.83 ± .40 


W+bb+Oj/Zbb+0] 
W+bb+lj/Zbb+1) 
W+bb+2j/Zbb+2j 


1.53 ± .03 
1.61 ± .04 
1.65 ± .06 



TABLE VII: Ratios of the cross sections for W + + N jets (including no b quarks) to W + + N jets (including two b quarks), 
and ratios of Z° /j* + N jets (including no b quarks) to Z° /j* + N jets (including two b quarks). Also given are the ratios 
W+ + bb + N jets to Z°/j* +bb + N jets . 
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FIG. 8: The ratios R for PDF X to R for CTEQ5L, where 
R = a(W + iVjets)/a(Z°/7* + Njets) . Changing the PDF 
affects this ratio much less - by at most 2% - than it affects 
the individual cross sections. 
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